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ABSTRACT 

Masking the horizontal branch and giant stars allows unambiguous measure- 
ments of mean age and metallicity in simple old stellar populations from metal 
and hydrogen line strengths. Billion year resolution is possible in the luminous 
halos of early type galaxies, constraining formation models. Most of the nui- 
sance parameters in stellar evolution are avoided by isolating the main sequence 
for analysis. The initial mass function and s-process element diagnostics may 
also be accessible. Integral field spectrographs have an significant advantage for 
this work, which is confusion limited by the presence of bright stars in medium 
to high surface brightness applications. 

Subject headings: galaxies: formation - techiques: spectroscopic - galaxies: 
abundances 



1. Introduction 

We know the age of the universe from the expansion rate and from the nuclear burning 
rate of the oldest stars. The agreement between these two approaches is one of the pillars of 
the standard cosmology (Mould 1999). To learn about galaxy formation we must be able to 
age date the halos of galaxies. This is difficult because a notorious degeneracy between age 
and metallicity affects the integrated spectra of galaxies (Faber 1973, Mould 1978, Schmidt, 
Bica & Dottori 1989). In this Letter we show that it is possible to deal with this problem by 
concentrating on the main sequence turnoff, where the age information in stellar populations 
resides, and by seeking separate spectral diagnostics of turnoff temperature and metallicity. 
Age resolution of a billion years is possible. Given the level of testing that main sequence 
evolution and the standard solar model have been subjected to, one can be reasonably 
confident about the uncertainties. The nuisance parameters of post main sequence stellar 
evolution: the mixing length to pressure scale height ratio, the Reimers mass loss parameter, 
the convective overshoot parameter and the horizontal branch (HB) color distribution are 
avoided. 
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A comprehensive study of the spectra of old stellar populations has been made by 
Worthey (1994). Mendel et al. 2007 point out that modelling of HB morphologies is difficult, 
as the interplay of contributing effects is not known well enough to be modeled purely from 
theory. Prescriptive methods have tended to be adopted. Maraston & Thomas (2000) 
find that a mass-loss prescription is able to reproduce the strong Balmer lines found in old 
elliptical galaxy populations, as well as the trends of increasing H/3 line strengths in low- 
metallicity Galactic globular clusters. Methods of breaking the HB/age degeneracy in the 
Balmer lines have been discussed by Schiavon et al. (2004). 



2. Isolating the main sequence 

Many of the galaxies we wish to study in a broad analysis of the formation of the first 
stars across the Hubble sequence of galaxies are at the distance of the Virgo cluster, where the 
main sequence turnoff of old stellar populations is at magnitude 35, intractable for current 
and planned telescopes. However, if we place our 30 meter class telescopes spectrograph slit 
across the outer parts of NGC 4472, say, and the seeing is good, we can in long exposures 
recognise individual horizontal branch stars at magnitude 32. If we study the light between 
these stars and further exclude stars with red colors (subgiants), we have a sample of turnoff 
stars, which are ideal for age and metallicity study. A quantitative example would be a study 
of a stellar population at 16 effective radii, 8.3 mag down from the central surface brightness 
of a de Vaucouleurs profile. If the local surface brightness is 27.5 mag/s 2 , a square arcsec 
would contain 2000 L Q at Virgo. This would be composed (by the fuel burning theorem, 
Greggio & Renzini 2011, Renzini & Buzzoni 1986) of some 500 main sequence stars and 4 
horizontal branch star^j]. This is a feasible proposition for a large telescope in natural seeing. 
Sky subtraction is a challenge at this surface brightness, of course. Working closer to the 
nucleus would require adaptive optics. 

An integral field unit spectrograph has a great advantage in this application. If the 
fibers are critically matched to the seeing, 4 fibers are lost in order to mask each horizontal 
branch star. However, in practise, rather than waste the light of HB stars, one would simply 
separate their spectra and make use of this information to constrain other parameters. The 
same goes for subgiant stars. 



*A referee has pointed out that a similar number of red giants would also be present and that analysis 
of their individual spectra would be valuable. Indeed this would be worth pursuing in order to measure 
detailed chemical composition. But it would not constrain age so well, which is our present emphasis. 
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3. Separating age and metallicity 

Consider two diagnostics, H/3 and the Ca K line, readily measurable in blue metal poor 
spectra. We write their line strengths, equivalent widths for example: 

P = P(g,T,Z) 

CaK = CaK(g,T,Z) (1) 

in the usual notation with g = log(gravity), T = log(effective temperature), and Z = 
log(metal fraction). Expanding around the main sequence turnoff we have 

6p= d^ 69+ &r 6T+ dz 6Z 

^ jz dCaK dCaK £rTi dCaK £ry 

5CaK = —6g + —5T + —5Z (2) 
We also have for turnoff stars 

r dg x . . dg 

69= m 6t + dz dZ 

dT dT 

ST =m St + az sz < 3 > 

where t is the log of the age of the simple stellar population. To evaluate these age and 
metallicity dependences, we use Padova YZVAR isochrones for helium fraction Y = 0.23 
(Bertelli et al. 2008). We find for Z = -4 that dg/dt = dT/dt = -0.0062 at the oldest 
ages and dg/dZ = 0.0945 & dT/dZ = -0.017. Synthetic line strength calculations from 
stellar atmospheres furnish the remaining information we need. According to Kurucz mod- 
elfiTigure 1) at [Fe/H] = -2, d(3/dT = 1.1, while at 6000K d(3/dg = -0.08 and d(3/dZ 
= 0.03. Similar results are obtained by Onehag et al (2008). We find from the synthetic 
spectrum calculations using MOOG (Sneden et al. 2012) and current ATLAS models (see 
Kurucz 1970, 1986 for the heritage): dCaK/dT = 50; dCaK/dg = 0.5 and dCaK/dZ = 
3. Ca K is measured as an equivalent width (EW) in A. MOOG calculates EWs by direct 
integration of the (continuum, line profile) difference. The EWs are already tabulated in 
the H/3 case. Based on the solar spectrum, one might question how well defined the CaK 
continuum is. However, there is no issue in the galaxy halo case Z-Z < -1. Solar abun- 
dance ratios are assumed and 2 km/s of microturbulence. Strictly one should also broaden 
the stellar line profiles with a galaxy velocity distribution. However, EWs are invariant to 
such post-broadening, e.g. for stellar rotation, instrumental broadening, or macroturbulence. 



2 http://wwwuser. oat. ts.astro.it/castelli/grids/gridni20k2odfnew/am20k2tab. html 
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Stellar population broadening has therefore not been done here. Similar results to Figure 2 
were obtained by Franchini et al. (2011). Thus we end up with 

5(3 = -0.08(-0.1245t + 0.09455Z) + l.l(-0.124<ft0.017<SZ) + 0.03SZ = -0.1265* + 0.0045Z 

SCaK = Q.lSt + 3.1SZ (4) 

One can formally eliminate metallicity or simply write: St ~ — 8Sf3. 

Now imagine comparing a 13 Gyr simple stellar population (main sequence only) with 
a 13 x 2.303St Gyr younger one. Subtracting the two spectra, one would see a S(3 given 
approximately by equation (4). Stars on the main sequence of lower mass than the turnoff 
would tend to cancel out. 



4. Secondary parameters 

In the preceding section the assumption has been that a stellar spectrum depends on 
(g, T, Z) and that age can be deduced from hydrogen and metal lines to the limits of 
observational error. However, the study of oj Cen has taught us that helium abundance is 
important and the systematics of [O/Fe] as a function of Z (e.g. Pilachowski & Armandroff 
1996, Cohen & Bell 1986) remind us that oxygen is the most abundant metal. Blue horizontal 
branch stars (Schiavon 2007) can tell us both helium abundance and oxygen abundance. In 
red giants the 6300A line measures oxygen abundance. Fibers that intercept the light from 
these stars can therefore be used to diagnose oxygen and helium anomalies in disintegrated 
spectra. However, this is beyond the scope of the present Letter. For the time being it is 
reasonable to assume that Z predicts O/Fe and that helium has its cosmological abundance 
in protogalactic halos. 

Alpha element enhancement has also been a subject of interest in the integrated light of 
early type galaxies. Figure 4 shows the prediction for the Mg/Fe ratio of metal poor turnoff 
stars. The iron line is 100-200 mA in equivalent width. Again, as a starting point in this 
Letter, it is assumed that ct/Fe is a constant in protogalactic halos and can be omitted from 
equation (1). Figure 4 is a precursor of a more detailed study, which, fortunately, is made 
easier by the existence of appropriate Kurucz models. The work of Dotter et al. (2007) 
notwithstanding, if a/Fe can be omitted from equation (1), it can also be omitted from the 
age error budget, provided a non-a element metal line is also observed. 
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5. Summary 

An observational error of 0.005 in /3 corresponds to an uncertainty of 1 Gyr at 12 Gyr. 
This is useful for distinguishing chemo dynamical models of rapid initial collapse from those 
of assembly through mergers up to redshift 2. A full evolutionary synthesis model, following 
Leitherer et al. (1999) and Maraston et al. (2003), is now warranted to include light from the 
full main sequence, rather than just the turnoff itself. Other lines should also be considered, 
such as the Ca triplet and the Paschen lines in the near infrared. Foster et al. (2009, 2010, 
2011) have pioneered the study of the Ca triplet in halos and globular clusters. 

From the very beginning, studies of integrated light have sought to constrain the initial 
mass function (IMF) in early type galaxies (Spinrad & Taylor 1971). Disintegrating the 
light will strengthen the observational approach to this problem too. Indices such as the 
Wing-Ford band (Wing & Ford 1969, Wing, Cohen & Brault 1977, van Dokkum & Conroy 
2011) will be more sensitive to IMF, when only main sequence light is analysed. 

Barium lines are an indicator of s-process activity (secondary nucleosynthesis) in the 
earliest stars. Ba II 4554A has 77 mA equivalent width at (T e ff, g, [Fe/Hl) = (6000, 4, 
-2). The ratio of predicted equivalent widths at other effective temperatures^ is shown in 
Figure 3. With a 100 km/s velocity dispersion this will be a 5% depression in the continuum 
of a metal poor main sequence spectrum. High signal-to-noise disaggregated^ spectra will be 
required. Figures 3 and 4 show that s-process and a elements are accessible to the technique 
described here and the temperature dependence of these line ratios is modest in the region 
which contributes most of the light. 
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Fig. 1. — Dependence of H/3 on temperature, gravity, and metallicity. The reference point 
for the top plot is g = 4, [Fe/H] = -2, for the center plot is T e// = 6000K, [Fe/H] = -2, and 
for the bottom plot is T e// = 6000K, g = 4. 
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Fig. 2. — Dependence of Ca K equivalent width on temperature, gravity, and metallicity. 
The reference points are the same as for Figure 1. 
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Fig. 3. — Dependence of the ratio of Ba 4554 to Ca K equivalent width on temperature. The 
reference point is g = 4 and [Fe/H] = -2. This curve reflects the behavior of two chemically 
similar elements in the same ionization state and ground state excitation, separated on the 
curve of growth by 4 orders of magnitude in abundance. 




Fig. 4. — Dependence of the ratio of the equivalent width of the first line of the Mg triplet 
to Fe 5269. 54A on temperature. The reference point is the same as for Figure 3. The single 
point indicates an alpha element enhancement of 0.4 dex. 



